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Abstract 

The lepton-number violating process u^N — ► fi~fi + fi + X is studied for 
the first time in connection with Majorana neutrino masses of the second 
generation. The sensitivity for light and heavy Majorana neutrinos is investi- 
gated. The ratio with respect to the standard model charged current process 
is improved by some orders of magnitude if compared to previously discussed 
Majorana induced AL^ = 2 processes. Non-observation of this process in 
previous experiments allows to demand the effective mass to be (m^) 10 4 
GeV, being more stringent than previously discussed direct bounds, however 
still unnaturally high. Therefore, in the forseeable future, indirect bounds on 
effective masses other than (m ee ) will be more stringent. 
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1 Introduction 



Investigation of lepton-number violating processes is one of the most promising 
ways of probing physics beyond the standard model. A particular aspect of this 
topic is lepton-number violation in the neutrino sector, which in the case of massive 
neutrinos would allow a variety of new phenomena Q. This emerges immediately 
in case of Majorana masses of the neutrinos, which are predicted in most GUT- 
theories ||. For v e the searches for Majorana neutrinos mainly rely on neutrinoless 
double beta decay (Oi//3/3), resulting in an upper limit on the effective Majorana mass 
(m ee ) =| Y t U^ m m m 7]^' | of about 0.2 eV H, where llljyyi £1X6 the mass eigenvalues, 



77^ = ±1 the relative CP-phases and U em the mixing matrix elements. In general, 
there is a 3 x 3 matrix of effective Majorana masses, the elements being 

(m a p) = \(U diag(mi^f p ,m 2 ?7f !p ,m377f p )C/ T ) a/ 3| 

= \ J2 mmV^UamUfSm with a, (3 = e, fi, r. 



(1) 



In this paper we explore the possibility to learn about Majorana masses associated 
with the second generation. The process under study is muon lepton-number vi- 
olating (AL^ = 2) trimuon production in neutrino-nucleon scattering via charged 
current reactions (CC) 

z/ M iV -> W+fi+X. (2) 

The muonic analogy for the quantity measured in neutrinoless double beta decay 
reads (m^) =\ Y.U 2 m m m rj^' | and is investigated for both light and heavy Ma- 
jorana neutrinos. The relevant diagram is shown in Fig. [I], which also defines the 
kinematics. Alternative ways discussed in the literature to obtain direct information 
about (m^) are muon capture on nuclei |4j] and lepton number violating K-decays 
like K~ — > 7r + \T ji~ m [?J §. The experimental knowledge of effective Majorana 
masses other than the one measured in Qv(3(3 allows only rather poor limits. The 
best values obtained are from muon-positron conversion in sulfur (therefore sensitive 
to (m Me ) 2 ) and lepton-number violating i^-decays: 

a (32 S+ -^32 gr +) J 0.4 GeV (singlet) 

a( S2 S + fi- =^ P* + u tl ) ^ ^ \ 1.9 GeV (triplet) (3) 

r( YTK= =?dlP < 1 - 5 - 10 ' 4 (m w ) < 1.3 • 10 5 GeV. 

Here the experimental limits are taken from the PDG |§ and for the mass bounds 
the theoretical results given in [TtJ and |7J] are used (all ratios are proportional to 



(iTi^a) 2 ). For muon-positron conversion two results are given, depending on whether 
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the proton pairs in the final state are in a spin singlet or triplet state, respectively. To 
our knowledge, there are no direct limits on other elements of (m a p). Note that we 
are considering direct limits, i. e. measuring processes which are directly dependent 
on the respective quantity without making any further assumptions. We will show 
that using the scattering process (0) instead of the rare K decay allows to set a 
limit better by one order of magnitude on (m^), which is however still too high to 
give a physical mass matrix in Eq. ([!]). Indirect bounds, e. g. using unitarity of the 
mixing matrix and oscillation experiments will of course be much more stringent. 
Direct production of Majorana neutrinos heavier than 100 GeV has been studied 
for various collider types (e~e + , e/i, pp, pp, e~p) |LJ], [L2| with typical results of a 
few to some hundred events per year for high-energy and luminosity machines. 



2 Model and calculation 

Using the diagram shown in Fig. [I] plus its crossed version we get for the squared 
invariant amplitude three terms, each factorizing nicely in three parts. At the upper 
and lower vertex we have the standard V-A term. The contribution of the Majorana 
neutrino, i. e. the part W + W + — > /i + yU + is well known from the theory of Ovftfl. 
For the calculation one might follow the strategy in Kayser's textbook [jn| or use 



—7= U^ml^lal-U^mW 



(4) 



the Feynman rules from [13 1, we will do the former. From here on we refer to this 
part of the diagram as the like" process. We use for particles the standard 

Lagrangian 

g_ 

where 7± = |(1 ±75)- We denote the Majorana neutrino with v m , the muon with \i 
and Upm is an element of the unitary matrix connecting weak interaction eigenstates 
with mass eigenstates. For the O^/3/3-like contribution to the first matrix element 
we have: 



(5) 



Mi oc [i/ ro 7 p 7+C/* m // 

To bring this in a form suitable for inserting the fermion propagator we use the 
relation: 

(6) 



z/ m 7 p 7+/i 



Here [f means the charge conjugated spinor of the muon. For a given spinor ip 
charge conjugation has the properties: 



ifj c = Cip , 
C~\C = -7 



^0 = -^T C -1 

C- l l,C = 7 5 T = 75- 



(7) 
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In the standard Dirac notation C = 27270 is the charge conjugation matrix. Since 
i^jji is a Majorana particle it has the property = A^i/ m , A m being a phase factor 



in the field expansion of v m connected with the intrinsic CP parity r^j p , see e. g. 



fl3| . For the expansion in terms of spinors and creation/annihilation operators the 
following relations are valid: 



u m oc fu + X m f + v, /i oc fu + f v, 
jjL c oc fu+ f + v, JF oc f + u + fv. 



(8) 



Here / annihilates a particle and / creates an antiparticle. Using all the above 
equations the matrix element describing the like process can be written as 

(see Fig. [I] for the definition of the kinematics) 



Mi oc A^t/;^/i c 7 p 7_z/ m z/ m 7^7 + /i 



(9) 



q 2 -m r 



From here on we neglect the mass m m in the denominator. See below for the case 
when this is no longer allowed. The above is the matrix element one would have 
obtained for an intermediate Dirac neutrino and applying the usual Feynman rules 
with one outgoing /i + written with an u instead of a v (thus producing a scalar 
expression) and one 7+ replaced with a 7_. Assuming CP invariance, the term 
can be written as (see e. g. [[13], PH ) 



'Sri fim 



^2 KnU, 



*2 



Vm n CF U 2 



(10) 



thus defining the usual effective mass. The matrix element is therefore proportional 
to the effective Majorana mass, just as in 0v/3(3 and the other mentioned lepton- 
number violating processes. 

For the crossed diagram, described by I2 is replaced by g 2 and k 2 by k 3 . Finally, 
the interference term is given by 



M\M 2 oc v(fc 3 )7^7+7M M ( fc 2M&3)7a7-7/M&2 



'111 



which has a negative sign due to the interchange of two identical fermion lines. Using 
the identities v = —u T C~ l and u = Cv T this can be written in a form suitable for 
using the completeness relations for the spinors: 



v{h)^ u ^ + ^^u(k 2 ) = -u T (k 3 )C 1 ^ + j tl Cv T (k 2 ) 
= -u T (k 3 ) (-i^ + -i u ) T v T (k 2 ) = v{k 2 )i l a + iuu{h). 



(12) 



4 



Putting all the couplings and propagators together, the matrix element for scattering 
with a quark can be written as (using the Feymann-gauge for the W propagator) 



\M\ 2 = |.Mi| 2 + |.M2 
--\{m^)\HAG%M^ 



if 



+ 25R (m\M: 
1 1 



Mw 



^{V^-AVVtei + m M )} + m q )l a l-{h + m q ,)} 



+ 



-2 



1 

7 2 



i 

I 
i i 

<?2 <?2 



^{l^l-lvih - m (U )7 a 7 + 7 /3 (^ + m^)} 

Tr{7 M 7_7 I/ (^j - 771^)7^7+7/3 (fe + «ty)} 
Tr {7/37+7^(^3 + 771^)7^ 7_ 7i/ (^ 2 - m M )} 



(13) 



1 



1 



2 12 (pi -P2) 



1 1 

~2""^ 

<?2 92 



1 



(fci -k 2 )(k 3 • fc 4 ) + 

((fe 2 • fc 3 )(fci • fc 4 ) - (&2 • fetX^l • &3 



1 



(fc 2 • fe 4 )(fc 1 • fe 3 ) 

{ki ■ k 2 )(k 3 



k 4 )) 



m g and m q i are the masses of the scattered initial and final state partons, respec- 
tively. Coupling to an antiquark is identical to replacing fc 4 with p 2 . The two short 
traces describe the SM V-A vertices, the ones inside the square brackets are the 
like process. Averaging over the parton spin adds an additional factor 1/2. 
The long traces were computed with the MATHEMATICA jT6| package TRACER 
||17|| . As can be seen, in the low mass regime the matrix element is proportional 
to (m wt ) 2 . If we take a heavy Majorana neutrino into account, one has to include 
the mass in the propagator for q 2 and q 2 , which we neglected from Eq. (^|) on. A 
statistical factor of 1/2 due to two identical final state muons has to be included in 
order to avoid double counting in the phase space integration. 
We also performed the calculation for purely right-handed (RH) currents and ob- 
tained the exact same result with the exchange W — > Wr. As known, RH currents 
must occur — if they exist — strongly suppressed with respect to the left-handed 
ones. Since the W momenta are relatively small in comparison to Mw, the cross 
section is proportional to 

/ 9 \ 4 



a cx G' F M w ((ql - M w )(q 2 3 - M w ))' 



G, 



oc 
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M 2 W , 



(14) 



forcing the purely RH case to be some orders of magnitude under the purely left- 
handed case, since M Wr > 6M W ||. Here we assumed Ql = 9r = 9- 
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One could also consider a heavy right-handed Majorana neutrino as suggested by 



some left-right symmetric theories [18|, where leptons are arranged symmetrically 




in left-handed (LH) and RH doublets, i. e. 

" , (15) 

Here is a heavy Majorana neutrino in a weak leptonic current of the form 

jf = 7Z 7 a 7+iV M + /I7 a 7-^ + • • • + h. c. (16) 

where the dots denote non muonic contributions. We consider it in order to illustrate 
the general properties of process (Q) in a model independent way and to stress the 
fact that the greatest sensitivity is achieved for a Majorana mass of 1 to 10 GeV, 
independent of the exact form of the coupling to the W, see below. Furthermore it 
serves as a comparison to the results from |5|, ||, who also considered this possibility. 
In general all possibilities could contribute at the same time. Performing the same 
calculation as before we get for the A^-case in Eq. ( p~3| ) a replacement (7+ <-> 7_) for 
the trace describing the 0z//?/3-like process which leads in the end to a replacement 
(&4 <-> pi, k& <-> p 2 ). Again, antiquark scattering is obtained by replacing with 
p 2 in the quark amplitude. 

Also possible is the exchange of other hypothetical particles such as those from the 
plethora of SUSY. Anyway, if process (Q) could be detected, general arguments 
guarantee a Majorana mass term for the muon neutrino, just as the Schechter and 
Valle argument [19j does in the case for neutrinoless double beta decay for the 



electron neutrino. In pC§ this theorem has been generalized to supersymmetry 
demanding also a non-vanishing Majorana mass for the scalar superpartners of the 
SM neutrinos. The smallness of the cross section however makes a more detailled 
analysis in this case not worthwhile: One could in principle derive limits on the 
right-handed coupling and/or Wr mass but they would definitely not compete with 
bounds derived by other methods, in contrast to the bound we will derive on 
in the next section. 



3 Results and Discussion 

For the evaluation of the total and differential cross sections we wrote a Monte Carlo 
program calling the phase space routine GENBOD ||21|| . For the parton distributions 
we chose GRV 98 (MS) NLO @ at Q 2 = s = ( Pl +p 2 ) 2 = x 2 M 2 + 2xM p E V} where 



M p denotes the proton mass, E u the energy of the incoming neutrino and x the 
Bj0rken variable. We set Q 2 = Q min whenever Q 2 went under the minimal allowed 
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value of 0.8 GeV 2 . To get the averaged neutrino-nucleon cross section we assumed 
an isoscalar target and replaced up- and down quarks to get the parton distributions 
for the neutron. 

Before presenting the results we estimate the ratio with respect to the total neutrino- 
nucleon cross section. The typical suppression factor one encounters when dealing 
with Majorana instead of Dirac neutrinos is M/E in the matrix elements with M 
being the Majorana's mass and E its energy. For the ratio R of the cross sections 
we have therefore: 

R ^(,N^,-^X) |10-- forM = 170keV 

cr(z/ M N - /x-X) \E) f w \ 10" 5 for M = 5 GeV 

(17) 

where we took as a typical value E = 30 GeV. For heavy neutrinos the behavior 
changes significantly: Instead of M/E we have now M~ 2 (see below) and the ratio 
goes as 

R oc F ™ 2 p ~ 10~ 7 for M = 100 GeV (18) 

These ratios will of course be further suppressed by a very small phase space factor, 
which rises slightly with energy and turns out to be about 10 -7 , as well as by factors 
arising from bounds on the mixing with heavy neutrinos, see below. 
As expected, the cross section is tiny: If we use for the mass of the Majorana neu- 
trino the current limit from the direct muon neutrino mass search, m u =170 keV 



23| , we find that the cross section for energies in the range 5 . . . 500 GeV is of the 
order a(3/i) ~ 1 . . . 10 2 • 10~ 33 b, being 20 orders of magnitude lower than the total 
neutrino nucleon CC cross section of <7cc — 1 . . . 10 2 ■ 10 -14 b for the same energy 
range. The E v dependence of the cross section can be fitted as a quadratically 
polynom, i. e. cr(3/x, E v ) = a ■ E v + b ■ E 2 which has to be compared to the linear 
dependence of the total CC neutrino-nucleon cross section. If we assume that this 
behavior holds up to ultrahigh energies (which it does not due to propagator effects) 
the cross sections would be roughly equal for E v ~ 10 20 GeV, far beyond any rea- 
sonable scale. 

The scaling with (m^) 2 holds up to masses of about 1 to 10 GeV. The RH iV M 
produces a signal in the same order of magnitude. We plot the trimuon cross sec- 
tion in Fig. H] together with the total CC neutrino nucleon cross section of about 
0.7 • 10- 14 £„/GeV b, multiplied with 10" 20 . 

Despite the small values, the ratio of the trimuon process described here is signifi- 
cantly more sensitive on (m^) than other discussed processes: Abad et al. |7j get in 
a relativistic quark model for the decay K + — > 7i~fi + fi + a branching ratio of 2 • 10 -22 
while Missimer et al. M estimate the ratio of /i _ /i + -conversion via capture in 44 Ti 
with respect to a normal CC reaction to be 4 • 10~ 25 for a few hundred keV Majo- 
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m m = 7 GeV 


m m = 80 GeV 


E v 


R 


Rcor 


R 


Rcor 


25 


1.1 • io- 14 


4.4 • IO" 24 


1.1 • IO" 16 


4.1 ■ 10~ 21 


50 


6.2 ■ IO" 14 


2.5 ■ IO" 23 


8.5 ■ IO" 16 


3.1 ■ IO" 20 


100 


3.1 ■ IO" 13 


1.2- IO" 22 


6.3 ■ IO" 15 


2.3 ■ IO" 19 


250 


2.1 ■ IO" 12 


8.4 • 10~ 22 


8.7- IO" 14 


3.2 ■ IO" 18 


500 


7.4- 10~ n 


3.0 ■ IO" 20 


6.1 • IO" 13 


2.2 ■ IO" 17 



Table 1: Pure (R) and "f/^m-corrected" (R cor ) ratios of the process for masses of 7 
and 80 GeV and different neutrino beam energies in GeV. 



rana. Thus the process (0) is about two orders of magnitude closer to the relevant 
standard model process than previously discussed Majorana induced muon-number 
violating processes. 

Considering now the massive case, i. e. including the Majorana masses in the prop- 
agator, the squared amplitude will now be proportional to the sum 



(19) 



For the sake of simplicity we skip for the moment the factors rj^U^ and consider 
only one mass eigenvalue, which turns out to dominate the cross section when it 
has an appropriate value. First of all, the cross section as a function of mass will 
rise quadratically until the propagator takes over and forces a (mass) -2 behavior. 
This is displayed in Fig. ^| where we plot the total cross section for different neutrino 
energies. As can be seen the maximal value of the cross section as a function of 
mass is obtained in the range 1 ... 10 GeV, rising slightly with E v . The reason for 
that is that the integration over the neutrino propagator has its maximum in this 
range. This fact makes the greatest sensitivity independent of the coupling of the 
Majoranas to the leptons or Ws. One can show that the heavy right-handed iV M 
displays the same behavior as the left-handed Majorana case shown in Fig. |3], which 
underlines this fact. In Fig. f| we display the ratio with respect to the total CC 
neutrino-nucleon cross section for a mass of 5 GeV as a function of the incoming 
neutrino energy. Note that this light masses are ruled out [[| and that for higher 
masses (M ^ 100 GeV) the ratio scales with M~ 2 . 

Up to now all the numbers given were for U^ m = 1. In this case, a maximum 
of 7.4 • 10 -11 of the CC cross section would be reached for a Majorana with mass 
of about 7 GeV. A neutrino beam of 500 GeV, coming from a high energy and 
luminous /i + yU~-collider with 10 13 CC events per year could in principle produce a 
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few hundred of such events. 

However, there exist already strong constraints on the matrix elements U^ m from 
the data. The DELPHI collaboration |23 examined the mode Z — > Vu m and found a 
limit of \U^ m \ 2 < 2 • 1CT 5 for masses up to m m ~ 80 GeV. For larger masses analyses 
of neutrino-quark scattering and other processes yield \U^ m \ 2 < 0.0060 P5| . This 
pushes the best sensitivity range about a factor of 10 towards higher values. 
In Table [l] we show the ratios R with and without taking into account the limits 
given above for different energies and for Majorana masses of 7 and 80 GeV. As 
can be seen one cannot get closer than at most 10~ 17 , even for a 500 GeV neutrino 
beam. In [|12|] finite width effects were found to increase the cross sections for direct 
heavy Majorana production significantly. However, these effects show up for high 
center-of-mass energies and high masses so that in our kinematical and mass sensi- 
tivity region these effects should be negligible. 

Nevertheless, also in the massive case the improvement compared to existing num- 
bers is some orders of magnitude: Halprin et al. || find a (?7 Mm -corrected) BR 
smaller than 3 ■ 10 -27 for K + — > 7i^fi + fi + and S + — > 7i~fi + fi + for a universally cou- 
pled 5 GeV heavy neutrino and Ng and Kamal || get a ([/^-corrected) branching 
ratio of 1.3 ■ 10~ 25 for a 2 GeV right-handed Majorana coupling to the W as in Eq. 
(|T6|). This means, for E v = 100 (500) GeV and few GeV Majoranas, process (0) 
is up to 5 (7) orders of magnitude closer to the standard model CC process than 
previously discussed muon-number violating AL M = 2 processes, which are induced 
by Majoranas. Interestingly the highest BR for the mentioned K decay in || is also 
in the range of 1 to 10 GeV. 

Though the cross section is probably too small to detect this process in the near 
future, it still allows to set bounds on (m^) . Let us assume an upper limit on a 
process like (fj) of the order 10~ 5 of the standard CC process (otherwise it would 
have been observed already, see Section f|) and take an energy of E v = 100 GeV. 
Starting at small masses, i. e. a oc (m^) 2 , we find (m MjU ) ^ 10 4 GeV. This has to 
be compared to (m^) < 1.3 ■ 10 5 GeV as obtained from K-decays . 



What is now the significance of this (m^) bound? Obviously, in a three-neutrino 
framework all elements of (m a p) should be roughly in the same order of magnitude, 
therefore at most a few eV, not much higher than the limit for (m ee ). In scenarios 
with additional massive neutrinos one has to include other processes as restrictions 
for possible mass matrix models. The bounds for mixing matrix elements with 
heavy neutrinos are typically in the order of a few 10~ 3 to 10 -2 for z/ e , and v T) 
thus roughly the same for all three families. It is impossible to reconcile the bounds 
for (m ee ), (m eiJi ) and with these conditions. If one considers heavy neutrinos 

and allows (m ee ) to be as high as the other two one can in principle fulfill the other 
conditions, but one is lead to contradictions when taking flavor changing neutral 
current processes into account, see the Appendix. It is thus not possible to con- 
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struct a mass matrix (m a p) spanning 14 orders of magnitude from (m ee ) to (m^). 
Therefore, indirect bounds will be far more effective than direct ones. 

4 Experimental considerations 

Several experiments report the observation of trimuon events J27| , |28| , |29fl . The ob- 
served ratio of trimuon events (having a lepton number conserving ( h) signature) 

with respect to single charged current events is of the order 10 -5 . First thought to 
provide evidence for physics beyond the SM the explanation was soon given in terms 
of CC reactions with dimuon production via meson decay, radiative processes or di- 
rect muon pair production from subsequent hadronic interactions []30], |31|, |32|. An 
acceptance cut on muon momenta to be larger than about 5 GeV was applied by all 
experiments. 

To extract a ( — h +) signature several background processes in typical wide band 
neutrino beams have to be considered. Among them the most severe are lepton 
pair creation due to antineutrino contaminations of the beam (also having a ( — h 
+) signature) and charm production with an associated pion or kaon decay as well 
as overlaying events with beam muons. Furthermore, going to high momenta some 
misidentification in the charge might lead to additional background. 
The observables found in the past to be most suitable for distinguishing the men- 
tioned standard processes from new physics were the momenta of the muons, their 
two- and three-body invariant masses and the azimuthal angle distribution between 
the leading muon and the other two. The leading muon (/ii) was defined as the one 
which minimizes the sum of the transverse momenta of the remaining two with re- 
spect to the direction of the W, (W = V — jli). A complete listing of all relevant 
distributions is not our aim, however, for the sake of completeness, we plot the 
distributions of the muon momenta, which might be used to identify the process. 
From Figs. | to ^ it can be seen, that for the light (m MM )-case the two ji + have 
relatively low energy, while the \x~ from the V-A vertex has a broad spectrum with 
significantly higher energy. This is no longer valid for a heavy Majorana where the 
difference of the muon momenta is less clear, but is becoming larger with increasing 
neutrino energy. However, the like-sign muons have typically the same momentum 
distributions, which is an important experimental signature. It is a general feature 
that the momentum difference gets bigger when the energy E u is significantly higher 
than the mass of the intermediate Majorana. For mass and energy being equal the 
distributions are more or less identical. 

A similar search to the one described here could also be done with beams looking 
for the corresponding process z/^N — > ji~X, which though turn out to have a 
cross section lower by a factor of 9 to 10. 
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5 Summary and conclusion 



We investigated the reaction v^N — > fi~{i + {i + X at fixed target experiments medi- 
ated by light and heavy Majorana neutrinos. Using the fact, that no excess events 
were observed in past experiments at the level of 10~ 5 with respect to charged cur- 
rent events, we could deduce a limit of (m^) ^ 10 4 GeV. This is more stringent 
than other direct results discussed on this quantity, but obviously not reconcilable 
with other laboratory experiments. 

Some general properties of process (Q) were discussed: The largest sensitivity was 
found for heavy Majorana neutrinos in the region between 1 and 10 GeV because 
of the fixed target kinematics. This was pushed towards approximately 100 GeV 
due to existing limits on U^ m . This is relatively independent of incoming neutrino 
energy and independent on the precise form of the couplings, as can be shown with 
a right-handed Majorana. In general, process @ is closer to the standard model 
CC process by 2 (up to a few 100 keV mass) up to 7 (> 1 GeV mass) orders of 
magnitude than previously discussed Majorana induced AL M = 2 processes. We 
state again that our 10 4 GeV is the best direct limit, not the best achievable limit. 
One could consider various modifications of process (0) in order to constrain non- 
standard model parameters connected with the muonic sector. For the case of 0zy/?/5 

(') 

limits on some Yukawa couplings A^ fc , describing R-Parity violating SUSY effects 



were deduced |33| , the bounds being up to four orders of magnitude more stringent 
than the ones obtained from other processes. In addition, for muon capture in 44 Ti, 
extensions of the standard model were found [[| to have branching ratios some or- 
ders of magnitude higher than the Majorana case, so that it seems worthwhile to 
apply them to process (H) as well. 

The smallness of the cross section however makes such a detailled analysis not very 
worthwhile. This might change for the case of a neutrino factory from a muon stor- 
age ring with a large number of interactions. Work about this topic is in progress 
and will be presented in the near future. 

We concentrated our analysis on neutrino beams, especially v^. Since the beam 
energies are much higher than the lepton masses, the same arguments as described 
here would hold for other fixed-target experiments using charged lepton beams. 
However, new background processes have to be considered here. 
Furthermore, also a lepton-hadron collider such as HERA, which also has the ad- 
vantage of higher yfs can be used. The same strategy that lead to the bound on 
(m^fj) can of course be applied to infer quantities as (m^ T ) or (m TT ), for which no 
direct limits whatsoever exist. Taking the appropriate channels, no SM processes 
faking the signal exist. This has in the meantime been discussed in detail in [34]. 
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A Appendix 

The task is to construct a mass matrix with (m ee ) ~ 0, (m. e/1 ) ~ 2 GeV and (m^) 
10 4 GeV. For the matrix elements for mixing with heavy neutrinos holds p5| 



J2 \UeH\ 2 < 6.6 ■ KT 3 and J2 \ U »h\ 2 < 6-0 ■ 1(T 3 . (20) 

First, let us assume that U% H ~ and Ufa = 6 ■ 10 -3 . Then we find that tuh = 
1.7 ■ 10 6 GeV, leading to U e jj — 1-6 ■ 10 -5 which in turn leads to a contribution to 
(m ee ) of 4.1 • 10 5 eV, in contradiction to our assumption. In general we found no 
solution for the allowed parameters. 

Conversely, we might want to use the fact that there is a bound on heavy (m# > 1 
GeV) contributions from the Heidelberg-Moscow experiment of |35] 



E^ l ^<5-10- 5 TeV- 1 . (21) 

Ignoring the condition (m ee ) ~ allows to find parameters capable of obeying the 
(m efJ ) and (m^) limits as well as Eqs. (|20|) and fl2T|) . For example, U^ H = U^ H2 = 
KT 10 , U 2 eIl3 = 4 • 10- 13 and Ul Hl = 5 • 10" 3 , Ufa, = -5 • 10" 4 , U^ 3 = 10' 5 with 
mn x = 100 GeV, m H2 = 1 TeV and m H3 = 10 6 TeV. 

Then again we have FCNC processes like /x — > ej, which are sensitive on = 



Ufj,iU e imf. The experimental value of the branching ratio, BR < 1.2 • 10 11 
and the theoretical value from |37j gives m^7 < 1.23 GeV, which is not fulfilled by 
the choice given. 



References 

[1] K. Zuber, Phys. Rep. 305, 295 (1998). 

[2] For a review see: R. N. Mohapatra, Unification and Supersymmetry, 2nd edi- 
tion, Springer Verlag, 1992; P. Langacker Phys. Rep. 72, 185 (1981). 

[3] L. Baudis et al, Phys. Rev. Lett. 83, 411 (1999). 



12 



[4] J. H. Missimer, R. N. Mohapatra, N. C. Mukhopadhyay, Phys. Rev. D 50, 2067 
(1994). 

[5] A. Halprin, P. Minkowski, H. Primakoff, S. P. Rosen, Phys. Rev. D 13, 2567 
(1976). 

[6] J. N. Ng, A. N. Kamal, Phys. Rev. D 18, 3412 (1978). 

[7] J. Abaci, J. G. Esteve, A. F. Pachero, Phys. Rev. D 30, 1488 (1984). 

[8] L. S. Littenberg, R. E. Shrock, Phys. Rev. Lett. 68, 443 (1992). 

[9] Review of Particle Properties, C. Caso et al, Eur. Phys. J. C 3, 1 (1998). 

[10] M. Doi, T. Kotani, E. Takasugi, Prog. Theor. Phys. Suppl. 83, 1 (1985). 

[11] S. T. Petcov, Phys. Lett. B 139, 421 (1984); F. del Aguila, E. Laerman, P. 
Zerwas, Nucl. Phys. B 297, 1 (1988); W. Buchmiiller, C. Greub, Nucl. Phys. B 
363, 345 (1988); E. Ma, J. Pantaleone, Phys. Rev. D 40, 2172 (1989); J. Kogo, 
S. Y. Tsai, Prog. Theor. Phys. 86, 183 (1991); D. A. Dicus, D. D. Karatas, P. 
Roy, Phys. Rev. D 44, 2033 (1991); J. Maalampi, K. Mursula, R. Vuopionpera, 
Nucl. Phys. B 372, 23 (1992); J. Gluza, M. Zralek, Phys. Rev. D 48, 5093 
(1993); A. Datta, M. Guchait, D. P. Roy Phys. Rev. D 47, 961 (1993); A. 
Datta, A. Pilaftsis, Phys. Rev. D 50, 3195 (1994); A. Hoefer, L. M. Sehgal, 
Phys. Rev. D 54, 1944 (1996); J. Gluza, M. Zralek, Phys. Rev. D 55, 7030 
(1997); F. M. L. Almeida et al, Phys. Lett. B 400, 331 (1997); P. Panella, C. 
Carimalo, Y. N. Srivastava, hep-ph/9903253. 

[12] G. Cvetic, C. S. Kim, C. W. Kim, Phys. Rev. Lett. 82, 4761 (1999); G. Cvetic, 
C. S. Kim, Phys. Lett. B 461, 248 (1999), B 471, 471 (2000) (E). 

[13] B. Kayser, F. Gibrat-Debu, F. Perrier, The Physics of Massive Neutrinos, 
World Scientific, 1989. 

[14] H. E. Haber, G. L. Kane, Phys. Rev. 117, 75 (1985). 

[15] S. M. Bilenky, S. T. Petcov, Rev. Mod. Phys. 59, 671 (1987). 

[16] S. Wolfram, Mathematica, Addison-Wesley, 1991. 

[17] M. Jamin, M. E. Lautenbacher, Comp. Phys. Comm. 74, 265 (1993). 

[18] J. C. Pati, A. Salam, Phys. Rev. D 10, 275 (1974); R. N. Mohapatra, J. C. Pati 
Phys. Rev. D 11, 566, 2558 (1975); G. Senjanovic, R. N. Mohapatra, Phys. 
Rev. D 12, 1502 (1975). 



13 



[19] J. Schechter, J. W. F. Valle, Phys. Rev. D 25, 2951 (1982). 

[20] M. Hirsch, H. V. Klapdor-Kleingrothaus, S. G. Kovalenko, Phys. Lett. B 398, 
311 (1997), ibidB 403, 291 (1997). 

[21] F. James, CERN 68-15 (1968). 

[22] M. Gliick, E. Reya, A. Vogt, Eur. Phys. J. C 5, 461 (1998). 

[23] K. Assamagan et ai, Phys. Rev. D 53, 6065 (1996). 

[24] P. Abreu, Z. Phys. C 74, 57 (1997), C 75 580 (1997) (E). 

[25] E. Nardi, E. Roulet, D. Tommasini, Phys. Lett. B 344, 225 (1995). 

[26] H. Nishiura, K. Matsuda, T. Fukuyama, Mod. Phys. Lett. A 14, 433 (1999). 

[27] B. C. Barish et ai, Phys. Rev. Lett. 38, 577 (1977). 

[28] A. Benvenuti et ai, Phys. Rev. Lett. 38, 1110 (1977), Phys. Rev. Lett. 38, 1183 
(1977), Phys. Rev. Lett. 42, 1024 (1979). 

[29] M. Holder et ai, Phys. Lett. B 70, 393 (1977), T. Hansl et ai, Nucl. Phys. B 
142, 381 (1978). 

[30] V. Barger, T. Gottschalk, R. J. N. Phillips, Phys. Rev. D 17, 2284 (1978), 
Phys. Rev. D 18, 2308 (1978). 

[31] J. Smith, J. A. M. Vermasseren, Phys. Rev. D 17, 2288 (1978). 

[32] R. M. Barnett, L. N. Chang, N. Weiss, Phys. Rev. D 17, 2266 (1978). 

[33] R. N. Mohapatra, Phys. Rev. D 34, 3457 (1986); M. Hirsch, H. V. Klapdor- 
Kleingrothaus, S. Kovalenko, Phys. Rev. Lett. 75, 17 (1995), Phys. Lett. B 352, 
1 (1995); K. S. Babu, R. N. Mohapatra, Phys. Rev. Lett. 75, 2276 (1995); M. 
Hirsch, H. V. Klapdor-Kleingrothaus, S. Kovalenko, H. Pas, Phys. Rev. D 53, 
1329 (1996); H. V. Klapdor-Kleingrothaus, Proc. NEUTRINO 96, Helsinki, 
June 1996. 

[34] M. Flanz, W. Rodejohann, K. Zuber, Phys. Lett. B 473, 324 (2000); ibid. 480, 
418 (2000) (E). 

[35] G. Belanger et ai, Phys. Rev. D 53, 6292 (1996). 

[36] M. L. Brooks et al. (MEGA collaboration), Phys. Rev. Lett. 83, 1521 (1999). 



14 



[37] R. N. Mohapatra, P. B. Pal, Massive neutrinos in Physics and Astrophysics, 
World Scientific Lecture Notes in Physics, Vol. 41, World Scientific, Singapore, 
1991. 



15 




Figure 1: Feynman diagram for the considered process. It is q 2 = q\ — k 2 = 
Pi — ki — k 2 . For the crossed diagram fc 2 and k^ are exchanged and we denote the 
corresponding momentum of the Majorana neutrino with q 2 = q\ — k 3 = pi — ki — k 3 . 
For the W momenta holds: qi = Pi — k x and g 3 = k 4 — p 2 - 
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Figure 2: Total cross section for process (H) with a left-handed (solid) z/ m together 
with the total CC z/N cross section times 10~ 20 (dashed). The (effective) mass for 
the neutrino is (m^) = 170 keV. 
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Figure 3: Total cross section for a left-handed Majorana neutrino as a function of 
its mass for different neutrino beam energies. No limit on U^ m was applied. 
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Figure 4: Ratio for process (Q) with respect to the total CC z/N cross section for a 
left-handed Majorana of (effective) mass 5 GeV. No limit on U^ m was applied. 
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Figure 5: Differential cross section at E v = 25 GeV for the momenta of the three 
muons for the case of a left-handed Majorana with effective mass of 170 keV. k\ is 
the muon momentum from the standard V-A vertex for the incoming neutrino, k 2 
and /C3 are the muons from the 0^/3/3— like vertices. 
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Figure 6: Same as previous figure for a mass of 5 GeV and incoming neutrino energy 
of 100 GeV. No limit on U 2 was applied. 




Figure 7: Same as above for a mass of 80 GeV and incoming neutrino energy of 500 
GeV. Note that in this case the momenta of the two /i + are always larger than the 
momentum of the \x~ . No limit on U 2 was applied. 



19 



